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Abstract. We prove that arbitrary superharmonic functions and 
superparabolic functions related to the p-Laplace and the p-parabolic 
equations are locally obtained as limits of supersolutions with de- 
sired convergence properties of the corresponding Riesz measures. 
As an application we show that a family of uniformly bounded 
supersolutions to the p-parabolic equation contains a subsequence 
that converges to a supersolution. 



1. Introduction 

This work gives approximation results related to the stationary p- 
Laplace equation 

(1.1) A p u = div(|Vu| p ~ 2 Vu) = 
and the time dependent p-parabolic equation 

(1.2) d t u - ApU = 0. 

The results and arguments extend to a more general class of equations 
as well, but we restrict our attention to the prototype equations for 
simplicity. We have also deliberately decided to exclude the singu- 
lar case p < 2 from our exposition. The solutions of these equations 
form a similar basis for a nonlinear potential theory as the solutions 
of the Laplace and heat equations do in the classical theory. Equa- 
tion fll.2jl is also known as the evolutionary p-Laplace equation and 
the non-Newtonian filtration equation in the literature. For the elliptic 
regularity regularity theory we refer to [5], [9] and parabolic to [8]. 

In the nonlinear potential theory, so-called p-superharmonic and p- 
superparabolic functions are essential, see |9], [10J and |18| . They are 
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defined as lower semicontinuous functions obeying the comparison prin- 
ciple with respect to continuous solutions of the corresponding equa- 
tion. The main focus of this work is in the parabolic theory, but some 
observations may be of interest already in the elliptic case. In their 
definition, the superparabolic functions are not required to have any 
derivatives, and, consequently, it is not evident how to directly relate 
them to the equation. However, by |13j a £>-superparabolic function 
has spatial Sobolev derivatives with sharp local integrability bounds. 
See also [3J, |3] and Using this result we can show that every 

superparabolic function u satisfies the equation with measure data 

d t u - dw(\Vu\ p - 2 Vu) =n, 

where fi is the Riesz measure of u. A rather delicate, yet relevant, 
point here is that the spatial gradient of a superparabolic function is 
not locally integrable to the natural exponent p. Consequently, the 
Riesz measure does not belong to the dual of the natural parabolic 
Sobolev space. For example, Dirac's delta is the Riesz measure for the 
Barenblatt solution of the p-parabolic equation. 

Since the weak supersolutions belong to the natural Sobolev, they 
constitute a more tractable class of functions. If we have an increasing 
sequence of continuous supersolutions and the limit function is finite in 
dense subset, then the limit function is superparabolic. Moreover, if the 
limit function is bounded or belongs to the correct parabolic Sobolev 
space, then it is a supersolution. In fact, this gives a characterization 
of superparabolic functions. The approximating sequence of supersolu- 
tions is usually constructed through obstacle problems, see [9] and [TB] 
for the elliptic case and [13] for the parabolic counterpart. On the other 
hand, the standard approach in existence theory is to approximate the 
Riesz measure of a superparabolic function with smooth functions as 
in [3J, @], m and [H]. 

Our argument is a combination of these procedures. More precisely, 
we show that for a given superparabolic function, there is a sequence 
of supersolutions, such that the Riesz measures are smooth functions, 
which converges to the original function both in elliptic and parabolic 
cases. This also provides a useful tool for existence results, because 
our method gives a converging subsequence for any bounded family 
of supersolutions. This is needed, for example, the proof of Theorem 
4.3 in [13], and we provide a detailed argument here. In Appendix, 
we also provide a useful technical result related to the approximation 
of a measure in the dual of the Sobolev space. This result may be of 
independent interest. 
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2. The p-Laplace equation 

We will assume throughout the paper that p > 2. As usual, W 1,p (Vt) 
denotes the Sobolev space of functions in L P (Q), whose distributional 
gradient belongs to L P (Q). The space W l,p (Q) is equipped with the 
norm 

IMIw^cfi) = IMIi>(n) + II VuiiiP^). 

The Sobolev space with zero boundary values, denoted by Wq' p (Q), is 
a completion of C£°(f2) with respect to the norm of W 1,P (Q). 

Definition 2.1. Let Q be a bounded domain in MP. A function u G 
W^(Q) is a weak solution of (II. ip . if 

Vu\ p ~ 2 Vu ■ Vipdx = 

for all ip G Cq°(Q). Further, u is a weak supersolution if the integral is 
nonnegative for all nonnegative test functions ip. 

By elliptic regularity theory, weak solutions are locally Holder con- 
tinuous after a possible redefinition on a set of measure zero and weak 
supersolutions are lower semicontinuous with the same interpretation, 
see for example [9]. 

Our argument also applies to more general equations of the type 

divA{x, Vw) = 

with the structural conditions 

(1) x i — y A(x,£,) is measurable for all £ G M. n , 

(2) £ i — y A(x,£,) is continuous for almost all x G f2, 

(3) A(x, • £ > a\£\ p , a > 0, for almost all x G and all £ G R n , 

(4) |-4(x,0l < /3|£| p_1 for almost all x G and all £ G M n , 

(5) („4(a;,£) — A(x,i])) ■ (£ — rj) > for almost all x G f2 and all 
l,r] G M" with £ ^ 77. 

However, for simplicity, we only consider the prototype p-Laplace equa- 
tion in this work. 

Definition 2.2. A lower semicontinuous function u : Q — > (—00, 00] is 
p-superharmonic, if u is finite in a dense subset of Q, and the following 
comparison principle holds in every Q' <e Q: If h G C(Q') is a weak 
solution to ( II. ip in Q' d Q and w > h on <9f2', then w > /i in fi'. 

Clearly, min(w, v ) and aw + /3 are p-superharmonic if u and w are 
and a, (3 G M with a > 0. We refer to |9J and |18j for more on p- 
superharmonic functions. 
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There is a relation between weak supersolutions and p-superharmonic 
functions. Weak supersolutions satisfy the comparison principle and, 
roughly speaking, they are p-superharmonic, provided the issue about 
lower semicontinuity is properly handled: every supersolution has a 
lower semicontinuous representative. In particular, a lower semicontin- 
uous supersolution is p-superharmonic. However, not all p-superharmonic 
functions are weak supersolutions, as the example 



shows. This function is p-superharmonic in R n , but it fails to belong to 
W l( £ (R n ), when 1 < p < n. However, the function u is a weak solution 
in W 1 \ {0} and it is a solution of 

— A p u = 5, 

where the right-hand side is Dirac's delta at the origin. 

Next we recall the obstacle problem in the calculus of variations. 

Definition 2.3. For ijj G C(ft) H W 1,p (ft), consider the class J 7 ^ of 
functions v G C(fi) n W 1j> (Q) such that v - ip G W^ P (Q) and v > ip 
in Q. The problem is to find a function u G J 7 ^ such that 



Jn 

for all v G J 7 ^. The function if) is called the obstacle and u the solution 
to the obstacle problem in Q with the obstacle ip. 

By [9], there is a unique solution of the obstacle problem provided 
Fip is not empty, because J 7 ^ is convex. Continuity follows from the 
standard regularity theory, see [9j. 

For every p-superharmonic function u and Q' d Q, there is an in- 
creasing sequence U{, i — 1, 2, . . . , of continuous weak supersolutions in 
Q', converging to u. Such a sequence is found by means of the obstacle 
problem as follows: By semicontinuity, there is an increasing sequence 
ipi, i — 1, 2, . . ., of smooth functions converging to u pointwise. Let Uj 
be the solution to the obstacle problem with obstacle Then it can 
be shown that the sequence has the desired properties, see [PJ and 



An important consequence of this approximation is a local summa- 
bility result for p-superharmonic functions. This in turn implies that 
the solutions to the obstacle problem converge in a Sobolev space, and 
the existence of the Riesz measures of p-superharmonic functions. We 




if p 7^ n, 
if p = n, 




[18j. 
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record these facts in the following theorem. The example above shows 
that the bounds are sharp. 

Theorem 2.4. Let u be a p-superharmonic function in Q and let Ui, 

i = 1, 2, . . . , be a sequence of solutions of the obstacle problem as above. 
Then 

(1) u G L\ oc (Vl) and |Vtt| G L^ oc (Q) for every r and q such that 

n n 
1 < r < (p — 1) and 1 < q < (p — 1). 

n — p n — 1 

(2) There is a nonnegative Radon measure [i such that —A p u = fi. 

(3) The sequence u i; i = 1,2,..., converges to u in W X( £(VL) as 
i — > oo. 

The equation —A p u = fi means that 

/ \Vu\ p ' 2 Vu ■ Vtpdx = hpdfj, 
Jn Jn 

for all ip G C^°(Q). Here fi is the Riesz measure of u. Observe, that for 
p < n, the gradient of a p-superharmonic function u does not belong to 
Lf oc (p,), but however, it is a weak solution of the above measure data 
problem. If, in addition, u G W l( £ (fl), then the restriction of jj, belongs 
to W- 1 ' 13 ' ({!'), which is the dual of the space Wi*{tf), for every Q' d Q. 

For later use, we record the following lemma. The proof is based on 
a characterization of Sobolev spaces and their duals by means of Bessel 
potentials, and is given in Appendix IA1 

Lemma 2.5. Let fi G W~ l ' p ' (VL) be a positive measure. Then there is a 
sequence fi, i = 1, 2, . . . , of positive smooth functions such that fi — > \i 
in W~ 1,P '(Q) as i — >• oo. 

Now we are ready to prove an approximation lemma for weak super- 
solutions. 

Lemma 2.6. Let u be a continuous weak supersolution in Q, Q' <s Q an 
open set and let fi be the Riesz measure of u. Then there is a sequence 
Ui, i = 1,2, ... , of solutions to —A p Ui = fi, where the functions fi are 
nonnegative and smooth with the following properties: 

(1) Ui u in W 1,P {VL'), and pointwise almost everywhere in Q' , 

(2) Vui — > Vu pointwise almost everywhere in Q' and 

(3) /i fi in W~ l ' p ' {Vt') asi-^oo. 
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Proof. Since \x is the Riesz measure of u, we have 



l(/w}| 



ip d/x 



Vu| p ~ 2 Vu ■ Vipdx 



for every <p G C^(fi'). Because u G W /1 ' p (r2 / ), we conclude that fi G 
W~ 1,p (fi')- By Lemma T2. 51 choose a sequence (/i) of smooth, nonneg- 
ative functions converging to ft in W~ 1,p and let Ui, i — 1,2, ... , 
be the unique solutions of the boundary value problems 

-A p Ui = fi in Q', 
Ui-ue W^ p (fl'). 

Using u — Ui as a test function in the equations for u and Ui and 
subtracting the obtained equations, we have 

/ (\Vu\ p - 2 Vu - \Vu^~ 2 Vui) ■ V{u - dx = (fi - U u-u t ), 

where (•, •) denotes the duality pairing between W c ,' p (fi') and W" 1 *'^'). 
The left-hand side is positive, so we may take absolute values on 
the right-hand side. By an elementary inequality for vectors and the 
Sobolev inequality, we obtain 

|V(tt - Ui)\ p dx 

< I (|Vw| p ~ 2 Vw- \Vui\ p ' 2 Vui) ■ V{u-ui)dx 



(2.1) 

= I (A* — fh u ~ u i) I — II A 4 ~~ fiWw-^v'in') \\ U ~ U i\\wg' p (n') 

< c ll^ - fi\\w-^p'(n')\\^( u -Ui)\\ L p(n'), 
where c = c(n,p). From this, it follows that 

||V(u - Ui) II^J,/) < c\\n - fi\\ w -i,pi ( n>y 

This implies that itj — >• u in W /1 ' p (fi') as « — > oo. □ 

As a consequence of the previous result, we obtain a corresponding 
approximation result for p-superharmonic functions. 

Theorem 2.7. Let u be a p-superharmonic function in Q, fl ; <e Q an 

open set and let fi be the Riesz measure of u. Then there is a sequence 
Ui, i = 1, 2, . . . ; of solutions to —A p Ui = f] in Q' , where the functions 
fi are nonnegative and smooth with the following properties: 

(1) Ui u inW 1,q (Q!) with 1 <q < n(p — l)/(n—l), and pointwise 
almost everywhere in Q' , 
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(2) V«j — > Vu pointwise almost everywhere in Q' and 

(3) the functions fi converge weakly to fi in the sense of measures 
as i — y oo. 

Proof. Let Ui, i = 1,2,..., be continuous solutions to the obstacle 
problem converging to u in W 1,q (Q'). For each i = 1,2, . . ., we may 
apply Lemma [2.61 to find a solution of — A p Wj = fi in Q' such that 

ii ~ n 1 
\\Ui — Ui ym.p(n') < -. 

i 

Since 

\\U — Ui\\w^,v(Q') < \\u — Ui\\w^<i(a') + \\ u i — Ui\\\Yh<l(W), 

we see that (uj) converges to u in W 1,q (Q') as i — > oo. By passing to 
a subsequence, if necessary, we may assume that Vuj — > Vw pointwise 
almost everywhere in Q'. Recall that q < n(p — l)/(n — 1); thus the 
sequence (| Viti\ p ~ 2 Vui) is bounded in L r (Q') for 1 < r < n/(n — 1). In 
addition, the limit of the weakly convergent subsequence is |Vm| p_2 Vm 
by the pointwise convergence. This allows us to conclude that 



lim / (pfi dx = lim / |VMj| p 2 V«j -Vtpdx 
|Vu| p ~ 2 Vu- V^dx 
<p dfj, 

for every ip G Cq°(Q'), and the proof is complete. □ 



n' 



3. The ^-parabolic equation 

In order to discuss the parabolic case, some preparations are needed. 
Let Q be an open and bounded set in R n with n > 1. We denote 

n T = n x (o,t), 

where < T < oo. For an open set U in M n we write 

VL tlM = Qx (ti,t 2 ), 
where < ti < t 2 < oo. The parabolic boundary of fl^ t a is 

d p n tut2 = (on x [t 1: t 2 ]) u (H x {ti}). 

The parabolic Sobolev space Z/ p (0, T; W 1,p (f2)) consists of measurable 
functions w : Q T — > [— oo, oo] such that for almost every t G (0,T), the 
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function x i— > u(x,t) belongs to W 1,P {Q) with the norm 

(3.1) \W\\lp(o,t;WI,p(si)) = (^J {\u\ p + | Vw| p ) dx dtj < 00. 

A function u G L p (0, T; W lj> {n)) belongs to the space L p (0, T; W Q 1,P (Q)) 
if x i — ^ u(x,t) belongs to Wo' p (fi) for almost every t G (0, T). The lo- 
cal space £f oc (0, T; W^(Q)) consists of functions that belong to the 
parabolic Sobolev space in every Q' tl t <s Qt- 

Definition 3.1. A function it G Lf oc (0,T; W^{Q)) is a weak solution 
of (OD in f2 T , if 

(3.2) - / u^fdxdt+ [ \Vu\ p ~ 2 Vu-V(fdxdt = 

for all test functions ip G C^°(Qt)- The function it is a supersolution if 
the integral in (13.21) is nonnegative for nonnegative test functions. In a 
general open set V of IR n+1 , the above notions are to be understood in a 
local sense, i.e. u is a solution if it is a solution in every set Qt 2 ,t 2 <s V. 

By parabolic regularity theory, the solutions are locally Holder con- 
tinuous after a possible redefinition on a set of measure zero, see [8]. 
In general, the time derivative Ut does not exist as a function. This 
is a principal, well-recognized difficulty with the definition. Namely, 
in proving estimates, we usually need a test function tp that depends 
on the solution itself, for example tp = uQ where ( is a smooth cutoff 
function. Then we cannot avoid that the forbidden quantity u t shows 
up in the calculation of <p t . In most cases, we can easily overcome this 
difficulty by using an equivalent definition in terms of Steklov aver- 
ages, as on pages 18 and 25 of |S]. Alternatively, we can proceed using 
convolutions with smooth mollifiers. 

Definition 3.2. A lower semicontinuous function u : Qt ( — oo ; oo] 
is p-superparabolic in Qt, if u is finite in a dense subset of Qt, and 
the following parabolic comparison principle holds: If h G C(Q't lt t 2 ) is 
a solution of (jl.2p in Q' ti t ($= Qt and h < u on the parabolic boundary 
d p Q' tlM , then ft, < it in Q h ,t 2 - 

It follows immediately from the definition that, if it and v are p- 
superparabolic functions, so are their pointwise minimum min(it, v) 
and it + /3, G R. Observe that it + u and au, with a > 0, are not 
superparabolic in general. In addition, the class of superparabolic func- 
tions is closed with respect to to the increasing convergence, provided 
the limit function is finite in a dense subset. We refer to [10J for more 
information about p-superparabolic functions. 
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A lower semicontinuous representative of a weak supersolution is 
p-superparabolic, see [T7], but as in the elliptic case, not all p-super- 
parabolic functions are weak supersolutions. For example, consider the 



Barenblatt solution B v 



B p (x,t) 




p/(j>-l)\ (p-l)/(p-2) 

, t>0, 

t < 0, 



where A = n(p — 2) +p, p > 2, and the constant c is usually chosen so 
that 



/ B p {x, t) dx = 1 



for every t > 0. The Barenblatt solution is a weak solution of (II. 2p in 
the upper half space 

{(x,t) e R n+1 : x e R n , t > 0} 
and it is a very weak solution of the equation 

d t B p - A P B P = 5 

in M n+1 , where the right-hand side is Dirac's delta at the origin. In 
contrast with the heat kernel, which is strictly positive, the Barenblatt 
solution has a bounded support at a given instance t > 0. Hence the 
disturbances propagate with finite speed when p > 2. The Barenblatt 
solution is p-superparabolic in R n , but it is not a supersolution in 
an open set that contains the origin, because it does not belong to 
£foc(0,T;<?(fi)). 

The obstacle problem in the calculus of variations is a basic tool in 
the study of the p-superparabolic functions as in the elliptic case. 

Definition 3.3. Let ip 6 C°°(R n+1 ) and consider the class T^, of all 
functions v G C(0 T ) such that v G L p (0, T; W 1 > P (Q)), v = ip on the 
parabolic boundary of and v > if) in The problem is to find 
u G J 7 ^ such that 

/ / Vm| p ~ 2 Vm ■ (Vv - Vm) + (v- u)^-) dx dt 
Jo Jn ^ dt ' 

>- I \v (x, T) — u(x, T)\ 2 dx 
2 Jn 

for all smooth functions v in the class J 7 ^. In particular, u is a con- 
tinuous supersolution of ( II. 2p . Moreover, in the open set {u > if)} the 
function u is a solution of ( II. 2p . 
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For existence results for smooth sets Q, see |2], |6] and [IE]. Charac- 
terizations of the solutions of the parabolic obstacle problem have also 
been studied in [19|. 

As before, we may construct continuous supersolutions converging 
to a given p-superparabolic function in a space-time cylinder flt lt t 2 by 
employing the obstacle problem. For details, see |13| and |15j . Below, 
we will also need the fact that the solutions to the obstacle problem 
have a time derivative in LP (0, T; W~ l,p (O)), which is the dual space of 
L p (ti,t2] Wq' p (Q)). The latter fact is contained in the existence result 
proved in [6]. Here we use the fact that the obstacle functions are 
smooth. The following summability estimates are proved in [12j. See 
also [3J and 0]. 

Theorem 3.4. Let u be a p-superparabolic function in VLt, Q'a t 2 ^ 
be an open set and let Ui, i = 1,2, ... , be a sequence of solutions of the 
obstacle problem in Q' ti t2 as above. Then 

(1) u G L\ oc {VLt) and |Vw| G L 9 oc (Qt) for any r and q such that 
1 < r < p — 1 + p/n and 1 < q < p — 1 + l/(n + 1), 

(2) there is a positive measure fi such that 

d t u — A p u = n 

in Qt an d 

(3) the sequence Ui, i = 1,2,..., converges in L q (ti,t2;W 1,q (Q')) 
as i — > oo. 

The Barenblatt solution shows that these critical integrability expo- 
nents for a p-superparabolic function and its gradient are optimal. 

Theorem 3.5. Assume thatu is a p-superparabolic function in Qt, let 
Q! t t2 <s Qt be an open set and let fi be the Riesz measure of u. Then 
there is a sequence Ui, i = 1,2, ... , of solutions to 

d t Ui - A p Ui = ft 

in fl' t t , where the functions fi are nonnegative and smooth, with the 
following properties: 

(1) U{ — y u in L q (t l7 1 2 ; W 1 '*^')) with 1 < q < p - 1 + l/(n + 1) 
and pointwise almost everywhere in Q' tl t2 , 

(2) Vui — > Vu pointwise almost everywhere in Q! t t , 

(3) d t Ui G L p \t 1 ,t 2 ; W- l > p \VL')) for every i = 1,2, .' . . , and 

(4) fc^ fi in LP'itu^W-^'itt')) asi^oo. 

We use the following approximation result. As in the elliptic case, 
the proof is based on Bessel potentials, and is found in Appendix [A] 
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Lemma 3.6. Let fi G L p '(0,T; W' 1 *' {tt)) be a positive measure. Then 
there are smooth functions such that fi^-fiin LP (0,T; W~ 1,p (Q)) as 
i — > oo. 

With the previous lemma, we can prove the following convergence 
result for solutions of the obstacle problems. 



Lemma 3.7. The claim of Theorem \3.5\ holds when u is a solution to 
the obstacle problem with a smooth obstacle. In addition, the sequence 
converges in L p (ti,t 2 ] W l,p (Vt')). 

Proof. Let fi be the Riesz measure of u, Q' d Q and < ti < t 2 < T. 
Since the obstacle is smooth, we have 

and consequently 



1 0^)1 



ip d/i 



ud t ip dx dt 



Vu\ p ~ 2 Vu ■ Vipdx dt 



< 



+ llv«iii; { 1 niiit2) iiv^iu P( n' tlit2 ) 

< (ll9t w llip'(ti,t 2 ;W-i.p'(n')) + II^IIlp^^,^)) IMIz^i^w^n')) 
for every ip G C^°(Q' ti t ). This implies that 

with 

W^W Lp' (ti^W-^P 1 {Q.')) — ll^*' u lliy(ti,t2;W- 1 .P / (n')) + II^ W IIlp(^ t )• 

We choose nonnegative functions fi G C°°ift' t t ), i = 1,2,..., such 
that 

fi^fi in L^,^^ 1 ^')) 

as i — > oo. Let Ui, i — 1, 2, . . ., be the unique solutions of the boundary 
value problems 

d t Ui - A p Ui = fi in 
Ui = u on d p VL' tit , y 
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Using u — Ui as a test function in the equations for u and Ui, and 
subtracting the obtained equations, we have 

(d t (u - Ui),u - 

(3.3) +/ (\Vu\ p ' 2 Vu- \Vui\ p ~ 2 Vui) ■ (Vw- Vu { )dxdt 

Jn' 

= (A* - fi,u~ Ui). 
A formal computation gives 

(d t (u — , u — Ui) = I (u — Ui)d t (u — u^ dx dt 



n' f 



(u — Ui) 2 (x, t 2 ) dx > 0. 



A rigorous argument can be based on Steklov averages, see [8]. We 
conclude that 



(3.4) Jn' t 



*1>*2 



|V(w - Ui)\ p dxdt < (fi- fi,u- Ui) 

— c IIA* — fi\\LP'(ti,t2;W- 1 'P'(n'))\\ u ~ M «llLP(ti,i 2 ;U ; ' 1 ' !3 (n'))- 



Here we also applied an elementary inequality for vectors. By the 
Sobolev inequality we have 

||« - Wi||£p(n{ liia ) < c||V(u - Ui)\\ LPi n' titt2 ) 

with c = c(n,p). From (13. 4p and the Sobolev inequality above, we 
obtain 

||V(« - ^)|U P (^ iit2 ) < c||/i - fi\\^ { ;X,w-^'my 

The right-hand side tends to zero as i — > oo. Thus we obtain 

Ui^u in U , {t 1 ,t 2 ',W^ p {Q!)) 

as i — > oo. By passing to a subsequence, if necessary, we obtain the 
convergence almost everywhere. □ 

Proof of Theorem \3.5[ The proof is now similar to that of Theorem 12. 7\ 
let U{, i — 1, 2, . . . , be the solutions to the obstacle problem converging 
to u, and apply Lemma [3771 to pick the functions Ui, i = 1, 2, . . . , such 
that 

1 

— Wj||i,j>( tlit2 ;v^i,p(n')) — ~- 

The convergence of to w then follows by using the triangle in- 
equality. For the weak convergence of the functions fi, we note that 
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(| Vui\ p ~ 2 Vui) is bounded in L r (Q' t t ) for some r > 1, and by the point- 
wise convergence we have that the weak limit must be |Vm| p ~ 2 Vm; thus 
a computation similar to the one in the proof of Theorem 12.71 shows 
the weak convergence of the sequence (/j) to /i as i — > oo. □ 

4. A CONVERGENCE RESULT FOR WEAK SUPERSOLUTIONS 

As an application of the above results, we prove a pointwise conver- 
gence theorem for a bounded sequence of weak supersolutions of the 
p-parabolic equation in the slow diffusion case p > 2. This kind of ap- 
proximation is needed the proof of Theorem 4.3 in [14], and we provide 
a detailed argument here. The idea of the proof is to use Theorem 13.51 
to choose approximations of the supersolutions Ui, apply a compactness 
result to the approximations to obtain a limit function, and then show 
that also the original functions converge to the same limit. The advan- 
tage of this approach is that the Riesz measures of the approximations 
are known to be functions, so that norms of their time derivatives can 
be estimated on time slices. 

Theorem 4.1. Let Ui, i = 1,2,..., be weak supersolutions of (jl.2p 
in Qt such that \m\ < M < oo for every i = 1,2,.... Then, for 
a subsequence still denoted by {ui), there is a function u such that 
u,i — > u and Vui — > Vu pointwise almost everywhere in Qt o,s i —¥ oo, 
and, moreover, u is a weak supersolution of (jl.2p in Qt- 

Proof. Let fa be the Riesz measures associated to ttj, i = 1,2,... For 

&s 1: s 2 <s ®%i,t 3 ^ ^t, we have 

llVitilli^n^) < cM, and fa{Vt' ti t2 ) < cM 

by Caccioppoli's inequality, see the proof of Theorem 4.3 in [14J. By 
Theorem I3.5[ for each i, we may choose Vi such that 

d t Vi - A p Vi = fi 

where /j is smooth with the properties 

H/ilU 1 ^'' ) < 2cM, and - v i \\ L9 (t uh . w i,«(c i y < - 

% 

for some q with p — 1 < q < p. 

We aim at applying a compactness result of [20J to the sequence {vj). 
To this end, we need bounds for the gradients and the time derivatives 
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of the functions V{. We have 

IPilliaCti^W^'Cn')) <|| M i|U 9 (ii,t2;VK 1 >'J(C')) + \\ V i ~ u i lU^ti^W^Cfi')) 

< C (M+ || V« i || z ,(Q U )) + ^ 
<cM + l. 

Then we consider a bound for the time derivative. From the equation 
satisfied by Vi we see that 

(4.1) dtvt e L\Q': uS2 ) + L r ( Sl , s 2 ; W^iff)) 

where r = q/(p — 1) > 1, and functions in L 1 ^" ) are identified with 
distributions in the usual manner. We denote the functional in the 
second component in (14.1 j) by L and get 



\Dp\ 



Vvi\ p ~ 2 Vvi ■ Vtpdxdt 



>n 
< \\Vv. 



Thus L is a bounded functional in the dual space L r (si, S2] W~ 1,r (Q")), 
because of the L 9 -bound for Vvi. 

The estimation of the remaining term heavily uses the fact that /, 
is an L 1 -function instead of a measure, and can thus be estimated on 
almost every time slice. By the Sobolev inequality for s' > n, we have 



fi" 



fi(x,t)(p(x)dx <||/j(-,t)||Li(n»)||v?IU-(fi'') 

<c||/i(-^)||Ll(fi")ll V ^IL s '(fi") 

<c\\fi(-,t)\\ L l( n „)\\ip\\ wl<s l {n „y 

for all <p e Co°(^")- Thus 

\\fi(-,t)\\w-i>s(Q") < c\\fi(',t)\\ L i(a") 
for almost every t, S\ < t < S2, and 1 < s < n/(n — 1) since Cq°(Q) is 
dense in W ' s (^")- Integrating this estimate in time, we see that 

\\fi\\L^(s 1 ,s 2 ;W-^(n")) < c \\fi\\Li(n» liaa )- 
Combining these estimates, we arrive at 

||^i|Ui(s 1 ,s 2 ;TV- 1 ' Q! (fi")) < c ; 

where 

q n 



1 < a < min 



p — 1 ' n 
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By the estimates in the previous paragraph and Corollary 4 in [20J, 
(i>j) is compact in L 1 (f2' s ' i ). Thus there is a function u such that 
Vi — y u in ) for a subsequence. For the same subsequence we 

have 

\\u - Wi||xi(n» iS2 ) ^IIm-^IU 1 ^^) + IK-^IUm^W 

<||«-^IUi(n'/ 1>S2 ) + ?^o 

as z — >■ oo. Thus for a subsequence we get the pointwise almost every- 
where convergence of Ui to u as i — > oo. 

We may pass from convergence in Q" S2 d Qt to the full set 
by the usual exhaustion argument. An application of Theorem 5.3 in 
[15J shows that u is a supersolution, as well as the pointwise almost 
everywhere convergence of the gradients. □ 

Appendix A. Approximation in the dual of Sobolev space 

In this section, we prove Lemmas 12.51 and 13.61 We use the charac- 
terization of Sobolev spaces on W 1 by Bessel potentials. For functions 
u in the Schwartz class S of rapidly decreasing functions and a > 0, 
define the operator 

T a u = g a *u, 

where g a is the Bessel kernel with the property that its Fourier trans- 
form g a is 

= (i + ierr /2 . 

Thus T a has an inverse T_ a , with the Fourier multiplier 

(i + iei 2 r /2 . 

It is clear that both T a and T_ Q map S into itself, so we may define 
both operators on the class of tempered distributions S' in the usual 
manner. Note that the kernel g a is symmetric and consequently T a and 
T_ a are self-adjoint operators. 

Let 1 < p < oo. The Bessel potential space is now defined by 

a ' p {W n ) = {ueS' :u = T a f for / 6 L p (R n )}, 

and its dual is 

(R«) = { v e S' : v = T_ a g for g e L P '{W 1 )}, 
with the duality pairing 

(u,v) a>p - a , P > = / fgdx, 
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and the norms 

||w||a,p(Rn) = H/HipQKn) and |M|-a,p'(Rn) = ||fl , ||Lp'(TRn). 

For a = 0, 1, 2, . . ., we have 

a > p (R n ) = W a ' p (R n ) and - a ' p ' (W 1 ) = W~ a ' p> (R n ) 
with equivalent norms, see [1] and [21J. 

Proof of Lemma \2.b\ Let fi G W~ l,p ' (Vl) be a positive measure. We 
construct a sequence of positive, smooth functions f\ such that 

Hm - fi\\w-^p'{a) 

as z — )• oo. 

We extend all functions in Wq' p (Q) by zero to M. n , and we also extend 
H to W~ 1,p (M. n ) by the Hahn-Banach theorem. Hence it suffices to find 
smooth functions /j defined on R™ such that 

Hm — /illiy-^p'OR™) ~~ ^ ^ 

as z — )• oo. Through Bessel potentials, we find a function g G L p '(IR n ) 
such that 

/i = T_ x g. 

Define 

/i e = 7] £ * /I 

where rj e is the standard mollifier. Since mollification, being defined by 
convolution, is multiplication on the Fourier transform side, and the 
same is true for T_ 1; we have 

(A.l) fJ, E = r) e * fx = rj £ * (T^g) = T^r), * g). 

Thus 

M — Me II VK-^p'fR") — C llM ~~ AM|-i,p'(r™) 

=c||r_i(fli -rj E * g)\\-i, P > {Rn) 

=c\\d ~Vs* 9\\lp'(R") -> 

as £ — >■ 0. Finally, we note that since // is a positive functional, the 
functions /i £ are positive, by the properties of the standard mollifiers. 

□ 

Proof of Lemma \3.b\ For the parabolic case, let us denote 

V = L P (R; W^iW 1 )) and V' = L P '(R; W~ 1,p ' (R n )). 

By the same extensions as before, it suffices to consider these spaces. 
Through Bessel potentials, we see that for \x G V, we have 

fi(t) = T-ig t (x) 
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for almost all t G R, where g t G L p '(IR n ). We define 

g(x,t) = g t (x), 
for (x,t) G M n+1 , and note that 

c llfi'llLP'(R n + 1 ) — IIHIv c ll5'llLP'(M n + 1 )- 

The approximation is then obtained by defining 

/J, £ = Ve* V, 

where the standard mollification is taken in IR n+1 . We use the Fourier 
transform in a fashion similar to (lA.ip . and get 

fjL e (t) = T-i[{rk*g)(t)] 

for almost all t, where T_i depends on the spatial variable only. From 
this we have 

HA* - fieWv < C \\9 ~Ve* flUxP^n+l) -> 

as e — > 0. For positivity, it once again suffices to note that the functions 
fi e are positive if the functional fi is. □ 
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